Landau Levels from the Bethe Ansatz Equations 
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The Bethe ansatz (BA) equations for the two-dimensional Bloch electrons in a uniform magnetic 
field are treated in the weak field limit. We have calculated energies near the lower boundary of the 
energy spectrum up to the first nontrivial order, ft corresponds to calculating a finite size correction 
for the excitation energies of the BA solvable lattice models and gives the Landau levels in the 
present problem. 
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The system of two-dimensional Bloch electrons in a 
magnetic field has very rich structure. The ratio of a 
lattice constant and a magnetic length, which are two 
fundamental length scales, is crucially important. The 
incommensurability brings fractal structures in the prob- 
lem which can be observed in the energy spectrum and 
the wavefunctions Famous Hofstadter's butterfly di- 
agram is a typical demonstration of the structure. It is 
fractal and has a self-similar structure. It also has an 
interesting relation to the quasiperiodic systems (qua- 
sicrystals) [Q. 

Recently intrinsic importance of a quantum- 
mechanical phase is stressed as a geometrical phase, 
which is essential in the Aharonov-Bohm effect, the 
Berry's phase, anyons and the quantum Hall effect 
The Bloch electrons in a magnetic field also have impor- 
tance as a stage of a typical realization of the geometrical 
phase. It was first noticed by Zak in the study of a mag- 
netic translation group |4| . The effect of the geometrical 
phase for the Bloch electrons was pursued in the studies 
of the quantum Hall effect and fundamental relations 
between the Hall conductance and several topological 
invariants were discovered (the Chern number |^ and 
the winding number on the complex energy surface [^). 
This is a geometrical aspect of the problem. 

For an algebraic property of the problem which origi- 
nates from the magnetic translation group, there was a 
breakthrough recently. Wiegmann and Zabrodin found a 
relation between the Hamiltonian of the Bloch electrons 
in a magnetic field and the quantum group Uq{sl2) 
Using the relation, the Schrodinger equation is rewrit- 
ten as a functional equation and the Bethe ansatz (BA) 
equations are derived. The quantum group is a rela- 
tively new mathematical concept. It is a kind of gen- 
eralization of the Lie algebra which is characterized by 
the so-called g-parameter. This g-parameter represents 
the non-commutativity of the elements in the quantum 
group which is related to the non-commutativity of the 
magnetic translation group. Mathematically rich struc- 
tures of the quantum group enable us to obtain new in- 
sights of the physical problem. Since the BA equations 
are high degree multi-variable ones, it is difficult to ob- 



tain an explicit solution. It usually occurs for the BA 
equations of the BA solvable lattice models (the XXZ 
chain and the Hubbard chain). However, for the present 
problem, the explicit solution of the BA equations for 
the zero energy was obtained by Hatsugai-Kohmoto-Wu 
[H,|o|. When the flux per plaquette is irrational, the 
distribution function of the BA roots is nowhere differ- 
entiable and the polynomial representation of the Uq{sl2) 
gives a polynomial with quasiperiodic coefficients. It im- 
plies that the wavefunction has a multifractal character. 
Although this explicit solution is only restricted to the 
zero energy, the other energies were also studied numeri- 
cally . Abanov, Talstra and Wiegmann gave beautiful 
results based on the string type ansatz when the flux per 
plaquette is an irrational golden mean Krasovsky 
derived an integral equation for the distribution function 
of the BA roots in the weak fleld limit and calculated the 
energy of the lowest energy band [Q. This corresponds 
to calculating a flnite size correction of the ground state 
energy in the usual BA solvable models. 

In this letter, we focus on the weak fleld limit 
where the flux per plaquette is given hy 4> — Bo^/^q in 
units of a flux quantum $0 — hc/\e\ (a is a lattice con- 
stant and i? is a magnetic fleld). Physically this weak 
field limit can be understood as a continuum limit from 
the lattice system where the lattice spacing a vanishes 
with the magnetic field B fixed. We calculate energies 
of the lowest several states analytically up to the lowest 
nontrivial order in cj). It corresponds to calculating finite 
size corrections for excited state energies of the usual BA 
solvable lattice models. Mathematically this weak field 
limit is a semiclassical limit where the classical SI2 is re- 
covered from the quantum Uq{sl2)- 

The Hamiltonian of electrons on a square lattice in a 
uniform magnetic field is given by 

H = T,+Ty + \i.c., 

where and Ty are the covariant translation 



operators — 

I]m,ncL,n-i-ie*^"'"Cm,«, and Cm,n IS an annihilation op- 
erator for an electron at a site [m,n). Discrete rotation 



1 



of the phases 0^ „ and 0,^ „ gives a flux per plaquette (j> : 

Spiaquotto ^m,n = ^m,n + ^m+l,n ~ ^m,n+l ~ ^m,n = " 27r0. 

In the diagonal gauge (0^,n — "'('^ + '^)'^) ^m,n — 
— 7r(TO+n+l)0), the Hamiltonian in the momentum space 
is given as 7J = Ek^(k), where H{k) = {e-'^-X + 
e^^y X^^)Y + Y-'^{e-'^y X + e'^- X-^) with 2Q x 2Q ma- 
trices X and r (Xij = 1 if i - j = 1 (mod 2Q), = 
otherwise, Y = diag(q, q^, . . . , q^'^ j). When one takes P 
to be odd and chooses momenta {k^, ky) on the so-called 
midband line = {kx + ky)/2 = tt/2 (mod tt/Q)), 
the Hamiltonian is given by 



= t{q-q-'){pc{B)+Pc{C)}, 



where Pc{B) and Pc{C) supplemented by pc{A), Pc{D) 
are the cyclic representations of the quantum group 
Uq{sl2) = {A,B,C,D} [QQ. Here the quantum group 
enters in the problem where the g-parameter is given by 
q = e"*^. Since we consider the case (f) = P/Q with mu- 
tually prime integers P and Q, it implies that the q is a 
root of unity {q^'^ = !)■ This cyclic representation is only 
possible for the q which is a root of unity. The classical 
sl2 is a spin algebra and it has a usual highest weight 
representation given by the differential operators. Corre- 
spondingly the quantum group Uq{sl2) also has a highest 
weight representation given by the difference operators. 
The highest weight representation implies that finite de- 
gree polynomials can span the bases. On the other hand, 
the cyclic representation does not have a classical corre- 
spondence. However, using a relation between the high- 
est weight representation and the cyclic representation, 
the discrete Schrodinger equations of the Bloch electron 
in a magnetic field are extended to a functional equation 
for a polynomial '^/{z) — Y\'k=i '^k) of the degree Q — 1 
M as 



i(z"^ + qz)'i'{qz) - i{z~^ + q~^ z)-^ {q~^ z) = e*(z). (1) 



Putting z = Zm and also comparing the coefficients of 
the highest order term, we have the BA type equations 
as 



qz^ + 1 



1 ' m=l,2 



fe=i 



q^k 



, ^, . . . , - 



e^iqQ{q-q-')Y,^k. (2) 



fc=l 



For the semiclassical case (p = 1/Q which we consider in 
the paper, all the BA roots are on a unit circle, \zk\ = 1 
(fc = l,2,...,g-l) (see FigJ^). 
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FIG. 1. Numerical solutions for the roots of the Bathe 
ansatz equations: (a) for the second lowest energy band and 
(b) for the third band (<;/> = 1/50). 

Therefore, we set Zk = e*"^" {tpk e R). Then the BA 
equations (^ are written as 



cos((p„ 



2Q' 



"5 1 sin|((^ 
M sini((^„ 



(3) 



In the limit Q —>■ oo, the distribution function p{(p) of 
the BA roots is well defined and the energy is given by 



2n 



e^^ p{uj)duj^ 



where p{(pk) = limg^oo 1/(54 and Sk = fk+i - <Pk- 
Recently, Krasovsky calculated the energy of the lowest 
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band up to the l/Q order and an integral equation for p 
was derived Extending the method in Ref. we 
calculate energies of the general bands in the vicinity of 
the lower boundary of the spectrum (e = —4 + 0(1/(5)). 
To obtain the energies, one needs the distribution of 
the BA roots. For the n-th lowest energy band with 
n = 0{1), Q — n roots are almost uniformly distributed 
on the left semicircle with a unit radius and the other n— 1 
roots are on the right semicircle jl^ (see Figjl]). We de- 
termine the arguments of the n—1 roots up to the order 
of unity. We denote the arguments of the Q — n roots 
and those of the other n—1 roots as 9 and r respectively 
{9 e (7r/2,37r/2) and r & (-7r/2, 7r/2)). In general, the 
n—1 roots consist of clusters separated by finite distances. 
The roots within each cluster converge to the single point 
in the Q cxd limit. Then the arguments of the BA roots 
are written as {rn, . . . , ti„i }, {r2i, . . . , T2„2}, • • ■ where 
Tii ci, T2j ^ C2, ■ ■ ■ [Q ^ oo) {a ^ Cj,i^ j). Further 
we assume c, ±f (see FigJI). Then the BA equation 
(0) for Til in the limit Q — )■ cxd gives 



= lim V In 



lim 

Q^oo 



E 

i=2 



In 



sm^ici-Ok - §) 
sini(Tii - Tu + §) 



sin ^(tii - Ti, 



(4) 



Note that if the first summation is naively estimated as 
integral, the integrand diverges at the points \m — k\ <^ Q 
[x <^ y means ^ ^ (Q ^ oo))- In order to avoid 
this, we divide the summation into two parts 5" and i?, 
where S is the contribution from k = m — N,...,m + N 
{1 <ti N Q) and R is that from the rest k (we take m 
such that 9m 7^ ±f {Q oo))- After canceling numer- 
ators and denominators in suitable pairs, S is estimated 

+ + "(q^) where gi < iV < Q 
are constants independent of Q. Since 



l+sin ■ 



we need another 



as — In 

and Vm and Prr 

we have two parameters ^ and ^ 
set of (similar) equations for the BA roots which can be 
derived from the functional equation (^ ||l^. From the 
two sets of the equations, we can fix the parameters as 
S — Pm — 1, s = 7r(p-|-i). Therefore the asymptotic 
distribution of the 9^ is determined as (see Fig.H) 



1\ TT 

2) Q 



Q 



{k^Q). 



The first term is estimated as In : , . Adding the 

l+sinci ^ 

BA equations for the roots {th, . . . ,ti„j^}, we have 



TLi In 



1 — sin ci 
1+sin ci 



0. It means ci = 0. Similarly, we have 



C2 = ■ • • = 0. Therefore, there is only one cluster on the 
right semicircle (see Fig.p and 
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(Q 



Now let us calculate the rest of the BA roots 9i and 
the energy for the {2p + l)-th lowest energy band (we 
assume that the Q is even for definiteness). By the sym- 
metry of the BA roots distribution P, p^ , p^ (see Fig.^ , 
we write arguments of the BA roots as = ^tTi ^fe = 
-7^ + f +E'=d ^-k - -Ok, A: - 1,2, . . . , § -p- 1 
where Aj = 9j+i ~ 9j — ^. Note that X]^=o vanishes 
in the Q ^ 00 limit, since Aj — o(^). Moreover, Afc 
is exponentially small for the k of the order Q {k ^ Q) 
@. Therefore, we set Ok ^ -n + ^ + + o(^) for 
k Q where s and S are constants independent of k and 
Q. Numerical solutions for 9k are shown in Fig.0. For 
the energy e, we have e — —4 + + ^ + "'(^' ^) from 
Eq.(H). To fix the energy, we calculate s and S. Now let 
us rewrite the BA equations (0) by 9^ and t; as 
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FIG. 2. Numerical solutions for the arguments of the 
roots 6k for the third lowest band [Q =50, 70, 100 and 200). 
The straight lines are the asymptotic behavior for Q 00, 
ek = -n + ik + 3/2)7r/Q (fc = 1, 2, . . . , (Q - 4)/2). 

Cases for the 2p-th bands are also discussed in the 
similar way. To summarize the results, the energy of the 
{n+ l)th lowest band (n = 0, 1, 2, . . .; n = 0(1)) is given 
as 



1\ 47r 
2, 



(5) 



In Fig.||, the five lowest energies are plotted as a func- 
tion of 1/Q (20 < Q < 200) and the straight lines are the 
asymptotic behavior in Eq.(||). 
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n = 0,l,2,...; n = 0{l), 
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FIG. 3. Numerical results for the energies of several 
bands near the lower boundary. Asymptotic behaviors are 
shown by the solid lines. 



We can also treat the general energy bands near the 
higher boundary of the spectrum. These (n+ l)-th high- 
est energies are given by ej^ = 4 — (n + + o(^) by 
the parallel argument. 

Now let us discuss the physical outcome of the results. 
In the absence of the magnetic field, one can recover the 
parabolic dispersion in the continuum model from the 
tight binding model as 



i?fc — —2t{coskxa + cos kytt) 
-U + ta^k^ = -At - 



2 1,2 



h'k 

2m* 



It implies that the effective mass is given by m* — 
n^l2aH. Since 1/Q = = BaV$o = \e\Ba^ jhc, we 
rewrite the above result (H) with this effective mass m* 
as 



2}? 1 



+ ?it^c 



(a^O), 



where lo^. = \e\B j m* c. These are usual energies of the 
Landau levels except the diverging energy shift. Now we 
can analytically recover the Landau levels from the BA 
equations of the Bloch electrons in a magnetic field . 
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